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Evidence for spin correlation in tt production 
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We present a measurement of the ratio of events with correlated t and t spins to the total number 
of tt events. This ratio / is evaluated using a matrix-element-based approach in 729 tt candidate 
events with a single lepton I (electron or muon) and at least four jets. The analyzed pp collisions 
data correspond to an integrated luminosity of 5.3 fb~^ and were collected with the DO detector at 
the Fermilab Tevatron collider operating at a center-of-mass energy ^/s = 1.96 TeV. Combining this 
result with a recent measurement of / in dileptonic final states, we find / in agreement with the 
standard model. In addition, the combination provides evidence for the presence of spin correlation 
in tt events with a significance of more than 3 standard deviations. 

PACS numbers: 14.65.Ha, 12.38.Qk, 13.85.Qk 



Although top and antitop quarks are produced unpo- 
larized at hadron colliders, their spin correlation can be 
measured, and a significant correlation is expected in the 
standard model (SM). The strength of spin correlation 
depends on the production mechanism and differs, for ex- 
ample, for qq and gg induced ti production [1]. Since the 
top quark decays through the electroweak interaction be- 
fore it can interact through the strong interaction ^ [3] , 
the spin orientation of the top quark at production is 
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reflected in the angular distributions of the final state 
particles f?^ . We present a measurement of the spin cor- 
relation of the t and t quarks to check its consistency 
with that expected in the SM. 

The tt spin correlation strength C is defined by 
d^crtt/(dcos6'idcos02) = cru{l — C cosOi cos02)/4, where 
att denotes the tt production cross section, and 9i, 02 
the angles between the spin-quantization axis and the 
direction of flight of the down-type fermion from the 
W boson decay in the respective parent i or t rest 
frame. It is related to the fractional difference A = 
{Na — No)l{Na + No) in the number of events Na where 
the top and antitop quark spins are aligned and those 
where the top quarks spins have opposite alignment, No, 
by C — A\aia2\ where is the spin analyzing power 
of the final state fermion under consideration. In next- 
to-leading-order quantum chromodynamics (NLO QCD) 
ai+ = 1 for the charged lepton in t £+i/^6 decays and 
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a J = 0.97 for the antidown quark in t — )■ dub decays [5]. 
The value A = +1 (—1) corresponds to fuUy parahel (an- 
tiparaUel) spins. Using the beam momentum vector as 
the quantization axis, the SM predicts Asm = 0.78tJ^-^^ 
at NLO QCD for pp coUisions at ^ = 1.96 TeV g]. 

Three tt spin correlation measurements based on the 
double differential angular distribution have been pub- 
lished so far [SHE]- However, none of them had sufficient 
sensitivity to distinguish between the hypothesis of spin 
correlation, as predicted by the SM, and no spin cor- 
relation. A fourth measurement was performed by the 
DO Collaboration in an analysis of 5.4 fb~^ of integrated 
luminosity in the tt dilepton channel, and reached an ex- 
pected sensitivity of 3 standard deviations (SDs) from the 
no-correlation hypothesis. In that analysis [9 , leading- 
order (LO) matrix elements (MEs) were used to measure 
the ratio / of events with correlated t and t spins to 
the total number of ti events by comparing Monte Carlo 
(MC) simulations with SM spin correlation and without 
spin correlation to data. 

In this Letter, we present the first measurement of the 
ratio / using the matrix-element approach in tt £-|-jets 
events. The t and t quarks are each assumed to decay 
into a W boson and a b quark, with one of the W bosons 
decaying directly or via a leptonic tau decay into an elec- 
tron or muon and the corresponding neutrinos and the 
other W boson decaying into two quarks. We use 5.3 fb^^ 
of integrated luminosity collected with the DO detector 
at the Fermilab Tevatron pp collider and combine our re- 
sults with the corresponding measurement in the dilepton 
channel [9j. 

A description of the DO detector can be found else- 
where [TU]. We use the same event selections as in the 
measurement of in the ^-fjets channel |11| . We re- 
quire one isolated electron with transverse momentum 
Pt > 20 GeV and pseudorapidity \ri\ < 1.1 [12|, or one 
isolated muon with px > 20 GeV and |?7| < 2.0, as well as 
an imbalance in transverse momentum > 20 (25) GeV 
for the e-|-jets (/i-fjets) channel. Events containing two 
isolated charged leptons with px > 15 GeV are rejected, 
to avoid overlap with the dilepton channel. In addi- 
tion, we require at least four jets reconstructed using 
a midpoint cone algorithm |13[ with radius TZ = 0.5, 
Pt > 20 GeV, and \ri\ < 2.5; the jet with largest trans- 
verse momentum must have pt > 40 GeV. Jets origi- 
nating from b quarks are identified using the output of a 
neural network where variables characterizing the prop- 
erties of secondary vertices and tracks with large impact 
parameters relative to the pp interaction vertex are com- 
bined [n]. 

The tt signal, with contributions from both qq — >■ ti 
and gg — >■ ti, is modeled using the mc@nlo [15 event 
generator with the CTEQ6M1 parton distribution func- 
tions (PDFs) |16|, assuming a top quark mass rrit — 
172.5 GeV. We generate ti MC samples both with and 
without the expected spin correlation, corresponding to 



A — 0.78 and A ~ 0, respectively \VQ. The events are 
further processed through HERWIG [TH] to simulate par- 
ton evolution, hadronization, and decays of short-lived 
particles, followed by a full detector simulation using 
GEANT P2! ■ We overlay events from random beam cross- 
ings on the MC events to model the effects of detector 
noise and additional pp interactions. The same recon- 
struction programs are used to process the data and the 
simulated events. 

The background can be split into two components: 
multijet background, where some of the products of 
hadronic partons are misreconstructed as an isolated lep- 
ton, and inherent background from SM processes with 
final states similar to that of the ti signal. In the 
e-|-jets channel, background from multijet production 
arises mainly when a jet with high electromagnetic con- 
tent mimics an electron. In the /i+jets channel it occurs 
primarily when a muon originates from the decay of a 
heavy-flavor quark {b or c) and appears to be isolated. 
The multijet background is estimated from data [llj . 
The SM background is predominantly from Ty-f-jets pro- 
duction, with smaller contributions arising from single 
top quark, diboson {WW, WZ and ZZ), and Z-t-jets 
{Z — !■ ee in e-|-jets or Z — > in /x-|-jets as well as 
Z tt) events. The W^-l-jets contribution is normal- 
ized to data using an iterative procedure, where the 
expected ti and smaller SM background contributions 
are subtracted from the data before application of 6-jct 
tagging [TT]. The differential distributions for ly-t-jets 
are taken from a simulation using the ALPGEN MC pro- 
gram |20j. All smaller SM background contributions are 
also estimated using MC simulations but normalized to 
their next-to-leading-order predictions. Diboson events 
are generated with pythia jH], single top quark pro- 
duction with the COMPHEP generator [52], and Z-|-jets 
events are simulated using ALPGEN. All MC background 
samples are generated using the CTEQ6L1 PDF 1161 . 
The evolution of partons and the hadronization process 
are simulated using pythia. A matching scheme is ap- 
plied to avoid double-counting of partonic event configu- 
rations [23]. 

To make optimal use of the kinematic information in 
ti events, we calculate signal probabilities Pggn for each 
event using the LO ME for the hypothesis of correlated 
{H = c) top quark spins, as predicted by the SM for 
qq — >■ ti, and for the hypothesis of uncorrelated {H = u) 
spins [H [21]. We can write Pggn as a function of the 
hypotheses H = c and H = u as: 

Psgn{x;H) =^^y /pdf(9i) /pDF(g2)dgidg2 
qiq2S 

with (Jobs being the LO qq ti production cross section 
including selection efficiency and acceptance effects, qi 
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and q2 denoting the fraction of the proton and antipro- 
ton momentum carried by the partons, /pdf represent- 
ing the parton distribution functions, s the square of the 
center-of-mass energy of the cohiding pp system, and d$g 
the infinitesimal volume element of the six-body phase 
space. Detector resolution effects are taken into account 
by introducing transfer functions W{x,y) that describe 
the probability of a partonic final state y to be measured 
as a; = (pi, . . . ,pn), where pi denote the measured four- 
momenta of the final state objects (leptons and jets). For 
the hypothesis H = c, we use the ME for the full process 
qq ^ tt — W~^bW~b — I'^vihqq'b [55| averaged over 
the color and spins of the initial partons, and summed 
over the final colors and spins [f . For the hypothesis 
H = u, we use the ME for the same process, neglect- 
ing the spin correlation. The total tt production cross 
section cr^j and the selection efficiency do not depend on 
spin correlation, thus the normalization factor dobs can 
be omitted in Eq. [l] To reduce the number of dimen- 
sions for the integrals, we assume the directions of the 
momenta of jets and charged leptons, and the electron 
energy are all well measured, and that the tt system has 
negligible transverse momentum. In addition, we use the 
known masses of the final state particles as constraints. 

As we use only four jets when calculating -Psgn, there 
are 24 possible jet-parton assignments. This further can 
be reduced to four when identifying the jets originating 
from b quarks. If more than two jets are ^-tagged, we 
select only the two jets with the largest 6-tag neural net- 
work probability as the b jets, and assume other jets to 
be light-flavor jets. Given the inability to distinguish the 
flavor of the two quarks from the W decay, as required 
for the definition of the spin correlation variable, both 
possible jet-parton assignments have to be considered in 
the Psgn calculation. Additional details of the Psgn cal- 
culation can be found in Ref . [26] . 

To distinguish between correlated and uncorrelated top 
quark spin hypotheses, we define, as in Ref. [9^, a discrim- 
inant R [17], 



R = 



Psg„(a;;i? = c) 



Psgn{x;H ^ c) 



(2) 



To measure the ratio /mcas of events with correlated spins 
to the total number of events, we form templates from 
distributions of R for tt MC events with and without spin 
correlation as well as background. Since the main sources 
of background are from multijet and W^+jets events, 
Psgn is only calculated for these two contributions. The 
smaller backgrounds are modeled using the templates for 
iy-|-jets production. The templates are compared to the 
distribution of R in the data, and the fraction /mcas is 
extracted through a binned maximum-likelihood fit. To 
minimize the dependence of the result on absolute nor- 
malization, we calculate the predicted number of events 
as a function of /mcas £^nd au, and extract both simul- 
taneously. Events used in the templates are required to 
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FIG. 1: The distribution of the discriminant R for ^+jets 
events with four jets and an invariant mass of the two hght- 
flavor jets within ±25 GeV of the mass of the W boson. The 
expectation (including background) for complete spin corre- 
lation as predicted by the SM (/ — 1) and the case of no spin 
correlation (/ = 0), as well as the ti contribution for /moas, 
where /moas was taken from the combined fit in the ^-|-jets 
and dilepton final states, are shown. The first and last bins 
include contributions from R < 0.37 and R > 0.60, respec- 
tively. The bin width is chosen to minimize the statistical 
uncertainty. 



have at least two fe-jet candidates; nonetheless, events 
with fewer than two 6-tagged jets are included in the fit 
to constrain the signal and background normalization. 
The fitting procedure and 6-jet identification criteria are 
the same as used in Ref. [28]. 

To enhance the sensitivity, we divide events into four 
subsamples as a correct jet-to-parton assignment greatly 
improves the discrimination power of R. The events are 
divided into two groups of events with exactly four jets 
and more than four jets to reduce the dilution from initial 
and final state radiation. To reduce the contamination 
from events in which a fo-quark jet is mistakenly taken to 
come from a W boson decay, these two groups are again 
separated according to whether the invariant mass of the 
two light-flavor jets is within or outside of ±25 GeV of 
the W boson mass. The ±25 GeV window is based on 
optimization through pseudoexperiments. The main sen- 
sitivity to spin correlation is obtained in the subsample 
with four jets and a dijet invariant mass close to the W 
boson mass, where the probability of selecting the correct 
jet combination is the highest. In Fig. [T] the measured 
discriminant R for the most sensitive sample is compared 
for data and templates of tt production with SM spin 
correlation and without spin correlation, including back- 
ground. 

We consider the same systematic uncertainties as used 
in the measurement of the tt production cross section [TT] 
and tt spin correlation in dilepton events '9]. These are 
included in the likelihood flt through free parameters. 
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where each independent source of systematic uncertainty 
is modeled as a Gaussian probabiUty density function 
with zero mean and an rms corresponding to one SD in 
the uncertainty on that parameter. Correlations among 
systematic uncertainties for different channels are taken 
into account by using a single parameter to represent the 
same source of uncertainty. 

We distinguish between systematic uncertainties that 
affect the yield of the signal or background and those 
that change the distribution of R. We consider the jet 
energy scale, &-jet energy scale, jet energy resolution, jet 
identification, 6-tagging efficiency and 6-jet misidentifi- 
cation rate, choice of PDF, and the choice of rrit in the 
calculation of Pggn as the uncertainties that affect the 
distribution of R. Systematic uncertainties on normal- 
izations include those on lepton identification, trigger re- 
quirements, the normalization of background, the lumi- 
nosity, MC modeling, and the determination of multijet 
background. We also include an uncertainty on the shape 
of the templates varying each template bin within its sta- 
tistical uncertainty. 

MC pseudoexperiments for different values of / are 
used to estimate the expected uncertainty on /moasj based 
on the maximum-likelihood fits that provide the depen- 
dence of / on /meas- The ordering principle for ra- 
tios of likelihoods ^\ is applied to the distributions of 
/ and /moas, without constraining /moas to the physi- 
cally allowed region. From a total of 729 events in the 




FIG. 2: (color online) Bands for 68%, 95% and 99.7% C.L. of 
/ as a function of /meas for the combined dilepton and ^-|-jets 
fit. The thin light-color line indicates the most probable value 
of / as a function of /moas- The vertical dotted black line 
shows the measured value of /meas ~ 0.85. 



TABLE 1: Summary of uncertainties on /meas for the com- 
bined fit in dilepton and ^-|-jets channels. 



Source +1SD -ISD 



Muon identification 





003 


-0 


003 


Electron identification and smearing 





009 


-0 


008 


PDF 





058 


-0 


051 


mt 





024 


-0 


040 


Triggers 





007 


-0 


008 


Opposite charge selection 





002 


-0 


002 


Jet energy scale 


A 

u 


fine; 


— U 


AOS 


Jet reconstruction and identification 





007 


-0 


035 


fe-tagging 





012 


-0 


012 


Normalization 





039 


-0 


043 


MC statistics 





015 


-0 


015 


Instrumental background 





003 


-0 


003 


Luminosity 





023 


-0 


023 


Multijet background 





007 


-0 


007 


Other 





007 


-0 


007 


MC statistics for template fits 





156 


-0 


156 


Total systematic uncertainty 





176 


-0 


184 


Statistical uncertainty 





251 


-0 


258 



^-|-jets channels with a ti signal purity of 90%, we ob- 
tain /meas = 1.15lQ;43(stat -I- syst) and can exclude val- 
ues of / < 0.420 at the 95% C.L. Since the samples of 
dilepton [5] and £+ijets final states are statistically inde- 
pendent, results from the two channels can be combined 
by adding the logarithms of the likelihood functions and 
repeating the maximum- likelihood fit. We obtain 

/meas = 0.85 ± 0.29(stat + syst) (3) 

and a ti production cross section of au = ^.IT^Qg^ pb, 
which is in good agreement with the SM prediction [3D] 
and previous measurements The statistical and sys- 
tematic uncertainties on /meas are given in Table [l] For 
an expected fraction of / = 1, we can exclude / < 0.481 
at the 95% C.L. For the observed value of /meas — 0.85, 
we can exclude / < 0.344 (0.052) at the 95(99.7)% C.L. 
We therefore obtain first evidence of SM spin correlation 
at 3.1 standard deviations. The probability to have a 
true value of / = for the observed value of /meas = 0.85 
is 0.16%. Figure [2] shows corresponding bands of con- 
fidence level. The ratio /meas can be used to obtain a 
measurement of the fractional difference Ameas by ap- 
plying it as a multiplicative factor to the NLO QCD 
prediction of Asm- ^meas = /meas X ^sM- This yields 
^mcas = 0.66 ± 0.23 (stat-Hsyst) [3T]- 

In conclusion, we have presented the first measurement 
of ti spin correlation using a matrix-element-based ap- 
proach in the £-t-jets channel. When combined with our 
previous result in the dilepton channel, we obtain signif- 
icant evidence for the presence of spin correlation in ti 
events with 3.1 standard deviations. 
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